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Abstract 

We use ab initio calculations to estimate formation energies of cation (transition metal) antisite 
defects at oxide interfaces and to understand the basic physical effects that drive or suppress the 
formation of these defects. Antisite defects are found to be favored in systems with substantial 
charge transfer across the interface, while Jahn-Teller distortions and itinerant ferromagnetism 
can prevent antisite defects and help stabilize atomically sharp interfaces. Our results enable 
identification of classes of systems that may be more and less snsceptible to the formation of 
antisite defects and motivate experimental studies and further theoretical calculations to elucidate 
the local structure and stability of oxide interface systems. 
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I. INTRODUCTION 


The remarkable electronic pr^erties of transition metal oxides, including high transition 
temperature superconductivity [l|, colossal magnetoresistance ^ and metal-insulator tran¬ 
sitions make them of fundamental importance for condensed matter physics. Interest has 
signihcantly increased following the fabrication of atomic-precision heterointerfaces which 


bring together different transition metal oxides with different bulk properties 


101. This 


ability to control materials at the atomic scale ho 
with entirely new properties and functionalities 


ds out the promise of creating systems 


16|. 


Realizing these exciting possibilities requires atomically precise interfaces. However, stud 


ies of interfaces separating simpler semiconducting materials 


17 


20l | show that antisite de¬ 


fects (exchange of atoms across the interface) may occur and can have crucial (and typically 
degrading) effects on near-interface electronic properties. In particular, emergent phenom¬ 


ena such as d-wave superconductivity and Weyl metal behavior are typica 
disorder, and clean samples are required for a convincing observation 21 


]y sensitive to 


231 1 . Although 


important work, in particular on defects at interfaces characterized by polar discontinuities. 


has appeared [2J-l27|, the subject of antisite defects at oxide interfaces has received relatively 


little attention. 

In this paper we consider antisite defects at AMO^/AM'O^ interfaces separating different 
members of the AMO^ class of pseudocubic perovskite transition metal oxides. In these 
materials the A-site is occupied by a lanthanide or an alkali earth ion (we consider A=La or 
A=Sr) and the M-site is occupied by a transition metal ion (we consider M, M' drawn from 
the first transition metal row). We focus on the situation in which the A-site is occupied by 
the same ion throughout and a change in the M-site ion defines the interface so an antisite 
defect corresponds to an exchange of M and M' ions across the interface 2^ . For all relevant 
combinations of M and M' we compute the defect formation energy, and then provide a 
physical understanding of the results in terms of the relative importance of charge transfer 
across the interface (leading to octahedral volume disproportionation that favors defects) 
and Jahn-Teller distortions (which inhibit defect formation). For metallic systems, itinerant 
ferromagnetism emerges as an additional factor inhibiting antisite defects. 

The rest of this manuscript is organized as follows. Section [TTl outlines the methods used; 
Section m presents our principal results, namely energies and local lattice structure for 
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different antisite defect combinations; Section |IV] gives an interpretation of the results in 
terms of charge transfer, structural distortions and itinerant ferromagentism; Section IVlis a 
summary and conclusion. 


II. COMPUTATIONAL DETAILS 


We perform density functional theory calculations 29|, l30| within the ab initio super- 


cell plane-wave approach M|, as implemented in the Vienna Ab-initio Simulation Package 
(VASP) 32|. We employ the Perdew, Burke and Ernzerhof (PBE) parameterization 33|] of 


the generalized gradient approximation (GGA) to the Kohn-Sham potential and projector 
augmented wave pseudopotentials [3J,l3^. The energy cutoff is 600 eV. We employ three dif¬ 
ferent types of simulation cells. For most of the calculations, we use a \/2 x \/2 x 2 supercell; 
to test the effects of inter-defect interactions we use a 2 x 2 x 2 supercell and to understand 
the effects of interface-interface separations a \/2 x \/2 x 8 supercell is used. A 8 x 8 x 6 
Monkhorst-Pack grid is used to sample the Brillouin zone of the \/2 x \/2 x 2 supercell. A 
5x5x5 Monkhorst-Pack grid is used to sample the Brillouin zone of the 2x2x2 supercell. 
A 8 X 8 X 2 Monkhorst-Pack grid is used to sample the Brillouin zone of the \/2 x \/2 x 8 
supercell. Both cell and internal coordinates are fully relaxed until each force component is 
smaller than 10 meV/A and the stress tensor is smaller than 10 kBar. Gonvergence of the 
key results was tested with a higher energy cutoff and a denser /c-point sampling and no 
significant changes were found. Gorrelation effects on the 3d orbitals are included using the 
VASP implementation of the rotationally invariant GGA+U approximation introduced in 
Ref. 36(]. We use U = 5.0 eV on the d orbitals for all the transition metal ions considered. 


For early transition metal ions (M=Ti, V and Cr), we choose J = 0.65 eV on the d orbitals 
and for late transition metal ions (M=Mn, Fe, Go and Ni), we choose J = 1 eV on the d 
orbitals. In order to shift the empty La 4/ states to higher energy, we also use Uls. = 9.0 eV 


on the / orbital, following the value used in previous work 


371. While the GGA-I-U method 


is only an approximate solution of the correlated electron problem posed by transition metal 
oxides, it is generally accepted as a robust method that captures the important trends in 
ground state energy and is computationally tractable, permitting surveys of wide ranges of 
interfaces. The most signihcant errors are in dynamical quantities that are not important 
for this work. 
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FIG. 1: Panel A: 2 x 2 x 2 supercell. Panel B: \f2 x \/2 x 8 supercell. Column 1: ideal interfaces 
with no antisite defects. Column 2: one antisite defect in the supercell. Column 3: two antisite 
defects in the supercell. Antisite defects are highlighted by the black ellipses. 

III. RESULTS 

A. Inter-defect and interface-interface interactions 

We define the defect formation energy as the difference between the energy of systems 
with and withont antisite defects. Compnting the energy of a single antisite defect at an 
isolated interface wonld require an inhnitely large computational cell. Practical calculations 
employ hnite supercells and therefore involve both a non-vanishing defect density and a hnite 
spacing between interfaces. To assess the degree to which our hnite supercell calculations 
are affected by non-vanishing defect densities, we studied a 2 x 2 x 2 supercell (40 atom in 
total, see Fig. [UAl) which can accommodate either one or two antisite defects (see Fig. [1] 
A2 and [T]A3). This corresponds to 25% or 50% defect concentration per interface. We 
restricted attention to ferromagnetic states to avoid issues of interplay between inter-defect 
spacing and magnetic ordering wave vectors. Results for three representative choices of A, 
M and M' are shown in panel A of Fig. O The symbols are the calculated energy differences 
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FIG. 2: Solid symbols (solid lines guide to the eye): total energies (with respect to the total energy 
of ideal interfaces) of various antisite defect configurations (see Fig. [T]) calculated A) in a 2 x 2 x 2 
supercell (40-atom) and B) in a \/2 x \/2 x 8 supercell (80-atom). Dashed lines: estimation of 
defect formation energies calculated in a \/2 x \/2 x 2 supercell (20-atom). N is the number of 
antisite defects in the supercell. 

between a system with N defects and a system with none; the slopes of the dashed lines give 
the formation energies estimated from the 50% defect concentration calculations. We see 
that using the higher defect concentration (50%) provides a formation energy which slightly 
overestimates that from the lower concentration (25%). The higher defect concentration 
(50%) can be accommodated in a smaller ^/2 x \/2 x 2 simulation cell. 

To assess the consequences of a hnite distance between interfaces, we study a \/2 x \/2 x 8 
supercell (80-atom in total) in which the two interfaces are separated by four unit cells. We 
consider three conhgurations: i) both interfaces are ideal (Fig. [Tpl); ii) one interface is ideal 
and the other interface has antisite defects (50% concentration per interface) and iii) both 
interfaces have antisite defects (50% concentration per interface). We compare in panel B 
of Fig. [2] results obtained on a \/2 x \/2 x 8 supercell to results obtained on a \/2 x \/2 x 2 
supercell. Isolating the interfaces does not change the sign of the energy difference but does 
somewhat increase the magnitude. These results indicate that a y/2 x \/2 x 2 supercell can 
be used as a conservative estimator of antisite defect formation energy. 
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FIG. 3: Sketch of atomic and magnetic structures considered in this work. Panels A: rocksalt 
configuration. Panels B: layered configuration. Large orange balls denote A-site ions La or Sr; 
small red balls denote O ions; shaded octahedra denote octahedra containing transition metal M 
(darker shade) and M' (lighter shade) ions, respectively. A1 and B1 are ferromagnetic ordering; 
A2 and B2 are checkboard G-type antiferromagnetic ordering; B3 are A-type antiferromagnetic 
ordering (spins are parallel in each layer and antiparallel between adjacent layers). Magnetic 
moments are schematically indicated by green/yellow arrows. 


B. Energetics 


We use the \/2 x \/2 x2 supercell to survey 21 LaMOa/LaM'Oa interfaces (M, M' = Ti,V, 
Cr, Mn, Fe, Co, Ni) and 15 SrMOs/SrM'Og interfaces (M, M' = Ti, V, Cr, Mn, Fe, Co). The 
\/2 X \/2 X 2 simulation cell is illustrated in Fig. |3l it consists of four perovskite primitive 
cells (20 atoms in total) and is large enough to accommodate both Jahn-Teller 38| and 


GdFeOa distortions 


as well as Neel antiferromagnetic ordering. The stacking direction 


of the layered structure is along the .2 axis. 

We estimate the defect formation energy for a given MM' combination as: 


■^formation — 


Emm' {L) 


( 1 ) 


where L is the conhguration of an ideal interface and R is the conhguration of one antisite 
defect, which for the computational unit cell used here corresponds to a rocksalt or double 
perovskite structure in which the M and M' ions populate alternate unit cells. We consider 
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FIG. 4: A: LaMOa/LaM'Oa. Al) DFT-calculated defect formation energy for each 

combination. The unit is eV per defect. A2) Structural parameters A and Q for each combination 
with the most favorable magnetic structure. A is shown in blue in the lower-right triangle and 
Q is shown in red in the upper-left triangle. The unit on the color bar for A and Q are %. 
A3) Comparison between the fitted defect formation energy and the DFT-calculated defect 

formation energy . AE^'^^ are obtained by minimizing D of Eq. (jH). B: SrMOa/SrM'Oa. 

Bl) Same as Al. B2) Same as A2. The red (blue) backslash denotes those combinations in which 
Smm/ is -1 (+!)• The other combinations without a backslash has a zero value of Smm/- B3) 
Same as A3. The solid squares denote that are obtained by minimizing of Eq. ([5|). The 

open squares denote AE^^^ that are obtained by minimizing D of Eq. (j3|). The combinations with 
a red or blue slash are explicitly labelled. 

both ferromagnetic (Fig. |2]A1 and Bl) and checkerboard G-type antiferromagnetic ordering 
(Fig.[2]A2 and B2). For the L conhguration, we also test A-type antiferromagnetic ordering 
(ferromagnetic planes with magnetization alternating between layers), since this magnetic 
ordering naturally hts the L conhguration (Fig. [2] B3). We always select the magnetic 
ordering that yields the lowest energy state. 

The energetics of antisite defects at oxide interfaces from our calculations are summarized 
in Fig.|3]Al and Bl (numerical results for energy differences together with information about 
the magnetic ordering and whether the system is metallic or insulating are provided in the 
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FIG. 5: Al) Top view of two vertically adjacent layers in the R configuration with a large oxygen 
octahedron (blue) and a small oxygen octahedron (purple). A2) Top view of two vertically adjacent 
layers in the L configuration with a naturally large oxygen octahedron (blue) and a naturally small 
oxygen octahedron (purple). Compatibility with the geometry of the L conhguration imposes 
strains (green arrows) on the system, compressing the large octahedra and expanding the small 
ones. Bl) Top view of two vertically adjacent layers in the L configuration with one anisotropic 
oxygen octahedron (purple) and one isotropic oxygen octahedron (blue). B2) Top view of two 
vertical adjacent layers in the R configuration with one naturally anisotropic oxygen octahedron 
(purple) and one naturally isotropic oxygen octahedron (blue). Compatibility with the geometry 
of the R configuration imposes strains (green arrows) reducing the bond disproportionation of the 
anisotropic material and inducing a disproportionation in the isotropic one. Rotations and tilts of 
oxygen octahedra are suppressed for clarity. 


Appendix). Blue indicates negative defect formation energies; for these cases we expect that 
the corresponding MM' hetero-interfaces are susceptible to antisite defects. Red indicates 
positive defect formation energies, suggesting those interfaces would be stable against defect 
formation. 














































IV. LOCAL LATTICE DISTORTIONS AND MAGNETISM 


A. Definitions 


In the previous section, we found that for many but not all M/M' combinations, antisite 
defects were favored. In order to gain insight into the factors favoring or disfavoring the 
appearance of antisite defects, we examine the correlation of the defect formation energy 
with other observables. 

A basic motif of the perovskite AMO 3 structure is the volume Vm of a MOq octahe¬ 
dron. Differences in octahedral volume between MOq and M'Oq octahedra are most easily 
accomodated in the R configuration, so we define the MOq, M'Oq volume difference as 


A 


MM' 


^ I Vm — Vm' I 

Vm + Vm> 


( 2 ) 


with the Vm and Vm' evaluated in the R configuration. 

A second important structural variable is a volume preserving Q 2 -type Jahn-Teller distor¬ 
tion of octahedron MOq in which one pair of M-0 bonds increases in length and the other 
pair decreases; both pairs of bonds alternate in the plane (see the left panel of Fig. [5]B1). As 
will be shown, the mismatch of Q 2 -type Jahn-Teller distortions (in-plane octahedral bond 
disproportionation) between MOq and M'Oq octahedra makes an important contribution 
to the stability of the L configuration. To quantify this, we define the bond-length dispro¬ 
portionation for ion M as Qm = 2|/i — hl/ih + h) where h and I 2 are the two in-plane M-0 
bond lengths for ion M in the L structure and then define the bond disproportionation 
mismatch Qmm' as 


Qmm' — \Qm — Qm'\ (3) 

for MOq and M'Oq in two adjacent layers in the L configuration. 

B. Analysis: La-based interfaces 

We begin our analysis with the LaMOs/LaM'Os interfaces, where for most MM' com¬ 
binations the ground state is insulating and the GGA+U method is expected to be reliable. 
Results for A are presented in a color scale in the lower right portion of the boxes in Fig.|4]A2. 
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Comparison to Fig. |4]A.l shows that a large octahedral volume difference is associated with 
a positive defect formation energy AE^^,. This correlation naturally arises from a strain 
effect. Fig. 041 shows that in the R configuration the large and small oxygen octahedra can 
be naturally accommodated. However, the geometry of the L configuration (Fig. 042) re¬ 
quires that the two oxygen octahedra have equivalent in-plane metal-oxygen bond lengths, 
inducing internal strain (represented by green arrows in Fig. 042) relative to the bond 
lengths preferred by the given charge configurations, thereby increasing the elastic energy 
of the L configuration. We note that although rotations of oxygen octahedra can accom¬ 
modate different octahedral volumes, our calculations on fully relaxed structures show that 
octahedral rotations can not reduce enough strain to favor the L configuration. 

For bulk LaMOs octahedral volumes change ^ 10% (the variation of lattice constant is 
^ 3%) as M is varied over the whole first transition metal row, but larger volume differences 
may occur in the superlattices. For example A = 15.1% for LaTiOs/LaNiOs [l3| and A = 
19.4% for LaTiOs/LaFeOs [^. These very large MOg volume differences are associated 
with complete charge transfer from M to M' ions, i.e. -|- M'^~^ —)■ (La- 

systems) or —)■ (Sr-systems). The MOg octahedral volume of 


electron acceptors expands while that of electron donors contracts. These substantia’ 
transfers are driven by large electronegativity differences between M and M' ions 


charge 


12 


m. 


A full list of the combinations in LaMOa/LaM'Os that have a complete charge transfer 
from M to M' is given in Table HI We see that interfaces at which significant charge transfer 
occurs are expected to be more susceptible to antisite defects. 

While charge transfer leads to octahedral volume changes that favor defects, the mis¬ 
match of Q 2 -type Jahn-Teller distortions between MOg and M'Og tend to inhibit defects. 
As can be seen from Fig. [5]B2, if a volume-preserving octahedral distortion has different 
amplitudes on the M and M' sites, it cannot naturally be accommodated in the R configu¬ 
ration (green arrows indicate strain). However, as Fig. [5]B1 shows, as long as the MO 2 and 
M '02 sheets have similar mean in-plane bond lengths, arbitrary layer-dependent Q 2 -type 
Jahn-Teller distortions can be accommodated without strain in the L configuration. Qmm' 
which mathematically defines the mismatch of Q 2 -type distortions in Eq. ([3]), is calculated 
for different MM' combinations using the L configuration and the results for Qmm' are 
presented in a color scale in the upper left portion of the boxes in Fig. |1]A2. Comparison 
to Fig. |1]41 shows that a large Jahn-Teller mismatch is associated with a negative de- 
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TABLE I: The combinations in La 2 MM '06 which have a complete charge transfer from M to M' 
ions (upper part) and which have a high-spin configuration (lower part). is the defect 

formation energy (Eq. (1) in the main text) using the most favorable magnetic ordering. The 
unit is eV per defect. Xmm' is the MOq octahedral volume difference using the R configuration 
(Eq. (2)). Qm describes the magnitude of Q 2 -type Jahn-Teller distortions of ion M using the L 
configuration (Eq. (3)). 


the combinations with a complete charge transfer 


charge configuration 

materials system 

A t^DFT 

Xmm' of the R configuration 

(E -(f - E 

La2TiV06 

-0.35 

23.9% 

E-E^E - E 

La2TiMn06 

-1.03 

27.0% 

E-E ^E-E 

La2TiFe06 

-1.24 

19.4% 

E-E ^E-E 

La2TiCo06 

-0.56 

18.4% 

00 

1 

t 

t- 

1 

As 

La2TiNi06 

-0.44 

15.1% 

E-E ^E-E 

La2VCo06 

-0.37 

20.4% 

t- 

1 

CO 

t 

to 

1 

■d* 

La2MnCo06 

-0.23 

24.0% 

00 

1 

CO 

t 

As 

1 

■d* 

"TS 

La2MnNi06 

-0.78 

18.0% 


the combinations with a high-spin configuration 


charge configuration 

materials system 


Qm of the L configuration 

E-E 

La 2 VMn 06 

0.26 

Qmh = 11.1%, Qv = 0.5% 

1 

CO 

La 2 CrMn 06 

0.20 

Qun = 10.3%, Qcr = 0.0% 

E-E 

La 2 MnFe 06 

0.28 

Qun = 14.9%, Qpe = 1.0% 


feet formation energy AE^^,. Among La-compounds, LaMnOs has the strongest Q 2 -type 
Jahn-Teller distortion j^. Our calculations confirm that the combinations LaVOs/LaMnOa, 
LaCrOa/LaMnOa and LaMnOa/LaFeOa, in which the Mn is in a high-spin state, all favor 
the L configuration, thereby tending to suppress antisite defects. We notice that in all the 
three cases, Qmii exceeds 10%, a value much larger than found for other ions (see Table [T]). 

As with the volume change, the relevant question for the Q 2 -typedistortion is the occu¬ 
pancy and spin state in a given structure, after any charge transfer has occurred. Usually 
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large ( 52 -type distortions are associated with negligible charge transfer, because charge trans¬ 
fer tends to create empty/half-filled/hlled d shells which are not Jahn-Teller active. Exam¬ 
ples include LaMnOs/LaNiOa {d'^ + d'^ —?■ d^ + d^) and LaTiOs/LaMnOa {d^ + d'^ —)■ d^ + d^). 
In both cases the charge transfer moves the Mn configuration away from the high-spin d^ 
state that favors Jahn-Teller distortions. 

The volume difference and Jahn-Teller effects will typically coexist and compete. To 
understand how this plays out in practice, we introduce a cost function 

^ — {aXuM' + /^Qmm')] (4) 

{MM') 

where the sum is over all the combinations MM'. Xmm' are calculated using the R con- 
hguration and Qmm' using the L configuration. Minimizing Eq. (|1]) yields a = —3.0 eV 
and = 2.7 eV for the La-based heterostructures. The opposite signs of a and (3 indicate 
that the volume change (AmmO and Jahn-Teller effect {Qmm') compete, as expected. The 
comparison between AE^^, and = aXMM' J- PQmm' is shown in Fig. |1]A3. While 

there is non-negligible scatter, the fit is reasonably good. In particular, the crude model 
correctly predicts the stability against defect formation (i.e. the sign of AE^^,) for most 
cases. 

C. Sr-based compounds 

Next we consider the Sr-based compounds. The defect formation energy AE^^jJ, for 
SrMOs/SrM'Os is shown in Fig. Hpl. The Xmm' and Qmm' for Sr-based heterostructures 
are calculated and displayed in Fig. |lp2. As for the La-based heterostructures, substantial 
charge transfer leads to large octahedral volume differences and favors the R configuration 
(examples include SrVOa/SrCrOs: dX—d? —)■ dX—dX and A = 5.6%; SrVOa/SrFeOa: d^—d‘^ —)■ 
dX — dX and A = 5.6%) while a large mismatch of (52-type Jahn-Teller distortions stabilizes 
the L configuration, for example SrFeOa/SrCoOe. However, we notice that the combinations 
(SrTiOa/SrFeOa, SrTiOs/SrCoOa, SrVOa/SrMnOa and SrCrOa/SrFeOa) that strongly favor 
the L configuration have a nearly vanishing Q 2 mismatch, suggesting the presence of an 
additional mechanism in the Sr-based heterostructures. 

We believe the additional mechanism acting in the Sr-based heterostructures is itinerant 
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ferromagnetism, which if present in the L configuration but not in the R configuration, sta¬ 
bilizes the L configuration. The MM' combinations with a ferromagnetic-metallic ground 
state in the L configuration and a non-ferromagnetic-metallic ground state in the R config¬ 
uration are labelled by a red slash in Fig. |1]B2. Comparison of Fig. |1]B2 to Fig. |3]B1 makes 
the stabilization effect evident. There is one case (SrCrOs/SrCoOs) where the R configura¬ 
tion is ferromagnetic metallic and the L configuration is not. As we will show below, in this 
case (labelled by a blue slash in Fig. |1]B2) itinerant ferromagnetism does not significantly 
contribute to the stabilization of the R configuration. 

To model the effects of metallic ferromagnetism, we define Smm' = 'S'mm' “ ^mm' where 
^MM' ^MM 'value 1 for ferromagnetic metallic states and 0 otherwise for the 
R and L configurations, respectively. We include this term in the cost function, obtaining 


fiSr — ^ ~ + I^QmM' + 7‘S'mM')] (5) 

(MM') 

Minimizing the cost flsr yields a = -1.8 eV, f3 = 2.0 eV and 7 = -0.4 eV. We comment that 
both in La-compounds and in Sr-compounds, a and [3 are very close in magnitude and of 
opposite signs, implying that the physical effects from the volume mismatch {Xmm') and the 
Jahn-Teller distortion mismatch {Qmm') are comparable and the competition between the 
two effects is a general phenomenon. The fit with Smm' is shown in Fig.|l]B3 as solid symbols 
and is clearly superior to the fit performed without Smm' (open symbols). Inclusion of Smm' 
makes the of the combinations with a red slash much closer to indicating 

a key role of itinerant ferromagnetism in stabilizing the L configuration. However, for the 
case with a blue slash, inclusion of Smm> drives further away from AE^^'^, implying 

that itinerant ferromagnetism does not significantly contribute to the stabilization of the R 
configuration. 

We make two comments concerning the Sr-based compounds. First, GGA-l-17 predicts 
charge ordering in a number of SrMOs/SrM'Os cases. Two-sublattice charge ordering is 
compatible with the \/2 x \/2 x 2 computational cell used here. It leads to an additional 
contribution to the energy difference that is not taken into account in Eq. ([5]). However, 
Eq. ([5]) works reasonably well, indicating that charge ordering does not substantially affect 
energetics. A more accurate many-body method (beyond GGA-I-G) is needed to study the 


delicate effects of long-range orderings 


41| . Second, SrFeOs, which in bulk is in the high- 
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spin df" configuration, is not Jahn-Teller active (unlike LaMnOa), presumably due to the 


enhanced covalency of Fe-0 bonding 


42|. 


V. CONCLUSION 


In conclusion, we have shown that M/M' antisite defects are energetically favored in 
wide classes of AMO^/AM'O^ heterostructures. The key driver of defect formation is a 
high degree of charge transfer across oxide interfaces, leading to large differences in equi¬ 
librium octahedral volume which in turn are most easily accomodated by antisite defect 
formation. On the other hand, a large mismatch of Q 2 -type Jahn-Teller distortion tends 
to inhibit antisite defects due to geometry constraints, as does itinerant ferromagnetism 
(in this calculation a signature of coherent metallic behavior across oxide interfaces). The 
association of defects with char ge t ransfer is unfortunate, as charge transfer is an important 
route to obtaining new physics 12l4l4|. 


Experimental 
reported 


43 


y, near-interface ‘dead layers’ of transition metal oxides are frequently 


44l |. suggesting a possible relevance of the present calculations to the be 


havior of real interfaces. On the other hand high quality LaTiOa/LaNiOs 1^ and 
LaTiOs/LaFeOs interfaces have been reported. Further experimental studies of transi¬ 
tion metal antisite defects at oxide interfaces would be very valuable for shedding light on 
this physics. On the theoretical side, we note that our calculations are based on the GGA-I-U 
approximation. While this method is believed to be a good approximation to the energetics 
of insulating systems, further investigation of selected cases using more sophisticated (but 
much more computationally expensive) methods such as dynamical mean held theory would 
also be desirable (although we emphasize that getting the local lattice structure correct is 
essential) and it is also interesting that our conclusions for the more metallic Sr compounds 
might also be revisited with other methods. We also note that we have used the same [/-value 
for all compounds. This choice is motivated by a desire to investigate chemical systematics 
without additional confounding factors but we note that our experience is that as long as 
U is not too small ^ 4 eV and not too large ^10 eV the basic physics of importance here 
(charge transfer, octahedral volume, Jahn-Teller distortions) are not particularly sensitive 
to U. 
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Appendix A: Energy difference, most favorable magnetic ordering and transport 
properties 

We show in Table UTI and Table HTTI the DFT-calculated energy differences AE^^,, defined 
in Eq. (1) in the main text. Table I is for La-compounds and Table II is for Sr-compounds. 
For each combination MM', we also show the most favorable magnetic ordering and trans¬ 
port properties for the rocksalt (i?) and layered (L) configurations. Text in blue (red) 
indicates those combinations that favor the rocksalt configuration (layered configuration). 
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TABLE II: La2MM'06 where M, M' = Ti, V, Cr, Mn, Fe, Co and Ni. The lower half is the energy 
difference between the R and L conhgurations AE = E{R) — E{L). The unit is eV per supercell 
(20 atoms). The upper half shows the most favorable magnetic ordering and transport properties 
for each conhguration. ‘F’ means ferromagnetic ordering, ‘G’ means G-type antiferromagnetic 
ordering, ‘A’ means A-type antiferromagnetic ordering. ‘M’ means metallic, ‘T means insulating. 
‘CO’ means charge ordering. 



LaTiOs (d^) 

LaVOs (d^) 

LaCrOs (d^) 

LaMnOs (d^) 

LaFeOs (d®) 

LaCoOs (d®) 

LaNiOs (d^) 

LaTiOg (rfi) 

bulk 

R: G-I 

L: A-I 

R: G-I 

L: G-I 

R: G-I 

L: F-M 

R: G-I 

L: G-I 

R: G-I 

L: G-I 

R: G-I 

L: G-I 

LaVOs (cP) 

-0.348 

bulk 

R: G-I 

L: G-I 

R: F-I 

L: F-I 

R: G-I 

L: G-I 

R: G-I 

L: G-I 

R: F-I 

L: G-I 

LaCrOg (d^) 

0.036 

-0.024 

bulk 

R: G-I 

L: A-I 

R: F-I 

L: G-I 

R: F-I 

L: G-I 

R: F-M 

L: F-M 

LaMnOs (d^) 

-1.028 

0.264 

0.196 

bulk 

R: G-I 

L: G-I 

R: G-I 

L: A-I-CO 

R: F-I 

L: F-M 

LaFeOa {d^) 

-1.244 

0.144 

0.144 

0.276 

bulk 

R: G-I 

L: G-I 

R: F-M 

L: F-M 

LaCoOs (d®) 

-0.556 

-0.368 

0.304 

-0.232 

-0.012 

bulk 

R: F-I 

L: G-I 

LaNiOs (d^) 

-0.444 

-0.264 

0.476 

-0.784 

-0.104 

-0.492 

bulk 
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TABLE III: St2MM'0q where M, M' = Ti, V, Cr, Mn, Fe and Co. The lower half is the energy 
difference between the R and L configurations AE = E{R) — E{L). The unit is eV per supercell 
(20 atoms). The upper half shows the most favorable magnetic ordering and transport properties 
for each conhguration. ‘F’ means ferromagnetic ordering, ‘G’ means G-type antiferromagnetic 
ordering, ‘A’ means A-type antiferromagnetic ordering. ‘M’ means metallic, ‘T means insulating. 
‘CO’ means charge ordering. 



SrTiOs {f) 

SrVOs {f) 

SrCrOs {f) 

SrMnOs {(f) 

SrFeOs {f) 

SrCoOs {f) 

SrTiOs {(f) 

bulk 

R: G(I) 

L: F(I) 

R: F(I) 

L: G(I) 

R: G(I) 

L: G(I) 

R: G(M) 

L: F(M) 

R: G(M) 

L: F(M) 

SrVOs {(f) 

0.056 

bulk 

R: G(I) 

L: F-CO(I) 

R: G(M) 

L: F-CO(M) 

R: G(I) 

L: A(M) 

R: F(I) 

L: G(M) 

SrCrOs {f) 

-0.072 

-0.644 

bulk 

R: G(M) 

L: G-CO(M) 

R: G-CO(M) 

L: F(M) 

R: F(M) 

L: A-CO(M) 

SrMnOs {f) 

-0.396 

0.476 

-0.088 

bulk 

R: F(M) 

L: F(M) 

R: F(M) 

L: F(M) 

SrFeOs {f) 

0.456 

-0.468 

0.244 

-0.184 

bulk 

R: F(M) 

L: F(M) 

SrCoOa {f) 

0.312 

-0.008 

-0.084 

-0.460 

0.268 

bulk 
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